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ABSTRACT 

Using microlensing measurements from a sample of 27 image-pairs of 19 lensed quasars we determine 
a maximum likelihood estimate for the accretion disk size of an average quasar of 4.0l3'^ light 

days at rest frame (A) — 1736A for microlenses with a mean mass of {M) = O.3M0. This value, in 
good agreement with previous results from smaller samples, is roughly a factor of 5 greater than the 
predictions of the standard thin disk model. The individual size estimates for the 19 quasars in our 
sample are also in excellent agreement with the results of the joint maximum likelihood analysis. 
Subject headings: accretion, accretion disks — gravitational lensing: micro — quasars: general 



^ . 1. INTRODUCTION 

C , The thin disk (Shakura & Sunyaev 1973) is the stan- 
dard model to describe the inner regions of quasars. This 
model predicts typical sizes for the accretion disks of 
: luminous quasars of ^10^^ cm (~0.4 light days) in the 

T-H ■ (observer frame) optical bands. However, recent mea- 
' surements based on quasar microlensing are challenging 

CN| this prediction. Based on a combined optical and X-ray 

^-H study of microlensing in 10 quadruply lensed quasars, 
Pooley et al. (2007) found that the source of optical 

. ^ light is larger than expected from basic thin accretion 
disk models by factors of 3 to 30. Studying microlensing 
\^ variability observed for 11 gravitationally lensed quasars 
, Morgan et al. (2010) also found that the microlensing es- 
timates of the disk size are larger (by a factor ^4) than 
would be expected from thin disk theory. Other mon- 
itoring based studies focused on individual objects like 
QSO2237-h0305 (Kochanek 2004, Anguita et al. 2008, 
Eigenbrod et al. 2008, Mosquera et al. 2009, Agol et 
al. 2009, Poindexter & Kochanek 2010), PG1115-I-080 
(Morgan et al. 2008), RXJ1131-1231 (Dai et al. 2010), 
HE1104-1805 (Poindexter et al. 2008) or HE0435-1223 
(Blackburne et al. 2011b) generally support these con- 
clussions. 

From single-epoch optical/IR and X-Ray measure- 
ments, Blackburne et al. (2()lla) used the wavelength 
dependence of microlensing (chromaticity) to study the 
structure of accretion disks in 12 4-image lensed quasars, 
finding disk sizes larger than predicted by nearly an or- 
der of magnitude. More detailed studies of individual 



objects based in microlensing chromaticity (Eigenbrod 
et al. 2008, Poindexter et al. 2008, Mosquera et al. 
2011, Mediavilla et al. 2011a, Muiioz et al. 2011) also 
found that the size estimated from microlensing is sub- 
stantially larger than the values predicted by the thin 
disk theory (by factors of about 10 in HE0435-1223, 5 
in SBS0909-I-532, and 4 in HE1104-1805). Upper lim- 
its on the sizes of MG0414-h0534 and SDSSJ0924-I-0219 
have also been set by Bate et al. (2008) and Floyd et al. 
(2009), respectively. 

Here, we will generalize these results by extending the 
study to the sample of 20 lensed quasars of Mediavilla et 
al. (2009, hereafter MED09). Our data consists of mi- 
crolensing amplitude measurements for 29 image-pairs 
from 20 lensed quasars based on a comparison between 
the flux ratios in the continuum and an adjacent emis- 
sion line. The emission line flux ratio provides an un- 
microlensed baseline that also removes the effects of ex- 
tinction (see e.g. Falco et al. 1999, Muiioz et al. 2004) 
and lens models in the determination of the magnifica- 
tion anomalies. As noted in previous works (Kochanek 
2004, Eigenbrod et al. 2008, Blackburne et al. 2011a), 
studies favouring the selection of objects with noticeable 
magnification anomalies (or epochs with high microlens- 
ing activity) may lead to a bias towards smaller quasar 
size determinations. The fact that one object does not 
show microlensing magnification anomalies may be due 
to the fact that it is either too big or that it lies in 
a region without significant magnification fiuctuations. 
Thus, non-detection of microlensing also puts constraints 
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on the size of the sources and should be taken into ac- 
count. The MED09 sample is also unaffected by this bias 
(indeed the histogram of estimated microlensing magni- 
fications in this sample peaks at Am = 0). However, 
the data from MED09 do not include, for most image- 
pairs, enough measurements at different wavelengths to 
address chromaticity. For this reason we will start dis- 
cussing the impact of chromaticity on the determination 
of microlensing based sizes (§2) to show that size esti- 
mates are dominated by microlensing amplitude and are 
relatively independent of chromaticity. In §3 we will use 
the microlensing data from MED09 to calculate a maxi- 
mum likelihood estimate of the accretion disk size. The 
estimate will be for the typical size of the accretion disk 
of an average quasar rather than that of any particular 
source. We also make individual size estimates for the 19 
objects using a complementary Bayesian approach, find- 
ing excellent agreement with the likelihood analysis. The 
main conclusions are presented in §4. 

2. SIZE DETERMINATION FROM SINGLE WAVELENGTH 
MICROLENSING MEASUREMENTS 

We will base our analysis on the microlensing mea- 
surements inferred by MED09 from optical spectroscopy 
available in the literature (29 quasar image pairs seen 
through 20 lens galaxies). The microlensing magnifi- 
cation between two images 1 and 2 is calculated as 

Am = (m2 - mi)micro = {m2 - mi)cont - {m2 - mi)line, 

where (m2 — rni)une is the flux ratio in an emission line 
and (m2 — rni)cont is the flux ratio of the continuum 
adjacent to the emission line. With this method the mi- 
crolensing magnification is isolated from extinction and 
from the mean lensing magnification, as these effects af- 
fect the line and the adjacent continuum equally. The 
consistency of this procedure has been confirmed with 
mid-IR data (MED09). 

The histogram of observed microlensing magnifications 
in MED09 showed that 93% of the quasar image pairs 
have microlensing |Am| < 0.8 mag. In this section, we 
are going to use this statistical constraint to estimate 
quasar accretion disk sizes from single wavelength mi- 
crolensing measurements. 

We will apply a statistical a pproach similar to that 
used bv lMediavilla et all ()2011aD to a representative lens 
system. Therefore, we take for the lens a standard SIS 
with values ki = 71 = 0.45 and K2 = 72 = 0.55 for the 
convergence and shear at the positions of the two images. 
These values are typical of the lenses in the MED09 (cf. 
their Table 4). For the redshifts of lens and source we 
have taken the values zi = 0.57 and Zs = 1.76, which are 
the average values for the sample. 

Using these parameters for the macro lens, we ge- 
nerate microlensing magnifications maps at the posi- 
tions of images 1 and 2 usi ng the inverse polygon 
ma pping algorithm d e scribed in iMediavilla et all ()2006f ) 
and iMediavilla et al.l H2011b[ ). The calculated maps are 
2000x2000 pixels in size, with a pixel size of 0.5 light 
days, or equivalently 1.295xl0^^cm. In the image plane, 
the covered area is the source plane area increased by fac- 
tors 1.5(1 — K ± 7)"^. For the fraction of projected mass 
density in stars, a, we use current estimates (Pooley et 
al. 2009, MED09, Mosquera et al. 2011) and consider 
only the cases with a = 0.05 and a = 0.1. We used 
microlenses of 1M0. 



In order to describe the structure of the source, we 
model the accretion disk with a Gaussian profile I{R) oc 
exp(— i?^/2rg), where the characteristic size is a wave- 
length dependent parameter (related to the half-light ra- 
dius as i?i/2 — l.lSrg). The wavelength dependence of 
the disk size is parametrized by a power law with expo- 
nent p such that rs{X) oc A^. We smooth the magnifica- 
tion maps with Gaussians scaled to two wavelengths As 
and Xn representing characteristic blue and red wave- 
lengths. We choose Xr/Xb = 2.3, the ratio between Lya 
and Mg II, which is roughly the wavelength range covered 
by observations in the optical band. 

We calculate the probability of observing microlensing 
magnifications Am = ^.2 — /^i (where /ii and /i2 are the 
magnifications at the positions of image 1 and 2 respec- 
tively) at wavelengths Xb and A^, namely, Amg and 
Amfl, for a grid in the size at Xb, Ts^Xb), and the ex- 
ponent p. We use a logarithmic grid in such that 
\ogr\ = 0.3 xi for i — - • • 11 {r\ in light days) and a 
linear grid in p such that pP = 0.25 xj for j = • • • 13. 
We restrict ourselves to the case |Atob| > |Ato/{|, as we 
assume that the size of the source increases with wave- 
length, and therefore the largest magnifications are ex- 
pected for the smallest (i.e. bluest) source^ 

The probability of observing Am^''* and Am^^^ in a 
model with parameters r\ and p^ is given by 

Prj,pj(AmB*, Am^^^) cx J dAniB J dAmRNije~^^\ 

(1) 

(2) 



where 

X 



2 _ (Ams - Am^^^)2 {Am,R - Am°^'f 



Here, Nij is the number of trials with Avib and Am^j 
for the case with parameters r\ and p' , and T^^obs and 
(J/^mobs are the uncertainties in the observed microlensing 
magnifications, for which we have taken a typical value 
of a^^^obs = cr^„o6= =0.05 mags. 

The calculations are performed for a grid of pairs 
AruBj^raR with Ams = 0.2xfc for k = —4- •■4 and 
Ato/j = 0.2 xm for m = 0---fc (so we always have 
\m\ < \k\ and \ AmB \ > \ AmR\) for a = 0.05 and a = 0.1. 
This range of values for the strength of microlensing 
covers the vast majority of observed cases as shown in 
the histogram of observed microlensing strengths in the 
sample of A4ED09 (their Fig. 1). Thus, the probabi- 
lity distribution p^'.pj (Am^''", Ato^***) is calculated for 
14x12x29x2 = 9744 cases. For each one of those cases 
(i.e. for every quadruplet (p, r^, Ato^*, Ato^^'*)), the 
probability is calculated using 10® trials. 

The results for the case a — 0.05 arc shown in Figure 
[1] for AniB > 0. We have shown only the right side 
of the full figure, as the cases with Aiub < are very 
similar to those with the same value of jAmsl so that the 
scenario is fairly symmetric with respect to Attib = 0. 
For the case a = 0.1 the results are very similar but 
the peaks in the probability distributions are displaced 
towards slightly higher values of r^. 

^ Inversions are possible but rare (see Poindexter et al. 2008) 
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Fig. 1. — Grid of two dimensional probability distributions p^i {Am'^" , Am'^") for the 15 cases with ArriB > and a = 0.05. Each 
column corresponds to the distribution for a given value of Amg indicated at the top of the column. Each row corresponds to the 
distribution for a given value of Amu indicated at the right. For each distribution the abcissac represent p from 0.0 to 3.25 and the 
ordinates represent logr^ from 0.0 to 3.3 (r^ in light days). The contour levels are drawn at intervals of 0.25o" (the contour at na is drawn 
at exp(— n^(T^/2) from the peak of the distribution). The contour at la is thicker than the rest. 



With the exception of the cases in which Am^ = 
Am/j, the distributions show a clear covariance in the 
sense that larger sizes imply lower values of the expo- 
nent p. We also see that larger values of Ams favour 
smaller values for the source size rs- 

The most important result from this figure is that the 
probability distribution with respect to the size is 
dominated by the strength of microlensing (Ams) with 
very little dependence on the amount of chromaticity. 
Conversely, the value of the chromaticity (Ato^ ~ ^'^ifl) 
is the main factor in determining the exponent p. In- 
deed, the most probable values of (values within la 
around the maximum of the distributions in Figure [1} 



remain roughly constant along each column in Figure 
[T] as we move through the different rows. That means 
that even if we lack information on chromaticity, we 
can still obtain valuable information on the size of the 
sources from measurements of microlensing at a single 
wavelength. Another interesting and robust result from 
Figure [T] is that, with the exception of the case with no 
detected microlensing {Ams = 0), the maximums of the 
distributions are all located in a fairly restricted range in 
Vs. Therefore, microlensing magnifications in the range 
< lAmsl < 0.8 will provide sizes in a rather restricted 
range between ~1a/M/Mq and ~16a/M/M0 light days 
irrespective of chromaticity, with larger magnifications 
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favouring smaller sizes. 

3. DISCUSSION 

Based on the results of the previous section, it is pos- 
sible to constrain the size of the sources from measure- 
ments of the microlensing strength at a single wave- 
length, even if we have no information on the amount of 
chromaticity present. As noted above (see also Kochanek 
2004, Eigenbrod et al. 2008, Blackburne et al. 2011a), 
each single microlensing measurement is affected by a 
degeneracy between the range of possible magnifications 
and source size effects. Statistics of large samples of 
lensed quasars (like MED09) can be used to minimize 
this uncertainty. In particular, we are going to use a 
maximum likelihood method to estimate an average size 
from the product of the individual likelihood functions 
for each image pair. This method includes the cases with 
little or no microlensing that by themselves would give 
rise only to lower limits on the size on an equal footing 
with those showing significant microlensing. 

The microlensing magnifications in this sample have 
been calculated using the continuum to line flux ratios 
of different lines for different objects. As size is related 
to wavelength, we should ideally use objects with magni- 
fications obtained from lines at similar wavelengths. On 
the other hand, to have good statistics, we should also 
try to keep the sample as large as possible. With these 
restrictions in mind, we have chosen a compromise in 
which we used all objects with magnifications measured 
in the wavelength range between Lya (1216A) to Mg II 
(2798A). With this choice, the average rest wavelength is 
in the relatively narrow range of (A) = 1736±373A, while 
we still keep 27 image pairs from 19 lensed quasars. The 
dispersion introduced by the effect of wavelength in the 
size estimate will be taken into account later. We then 
use the measured microlensing magnifications in these 
27 pairs as observed Attib (where B stands now for the 
average wavelength (A) = 1736 A). We computed magni- 
fication maps for each one of the images of the 27 pairs in- 
dividually taking the k and 7 values from MED09 (except 
for SBS0909+532 for which we take the values in Medi- 
avilla et al. 2011a). We again consider only a — 0.05 and 
a = 0.1. The calculated maps are 2000x2000 pixels in 
size, with a pixel size of 0.5 light-days. We have used the 
same logarithmic grid in Vs from §2, as well as a linear 
grid rl = 1.0 -|- 2.0xi hght days for i = • ■ • 24. Using 
the magnification maps, the probability of observing a 
microlensing magnification Arn°''* in image pair j for a 

model with parameter is given by: 



where 



2 {Am - Amf^f 



We can calculate a likelihood function for as 
i(r:)cx J]Pn(A™f^), 

where j runs over the 27 image-pairs considered. 



(3) 



(4) 



(5) 




log (r,/lt-day) 

Fig. 2. — Likelihood functions for the size of the accretion disk 
logrs for a = 0.05 (continuous line) and a = 0.1 (dashed line) in 
the full sample of quasars by MED09. 

The results of this procedure for a = 0.05 and a = 0.1 
are shown in Figure [21 The maximum likelihood occurs 
for logTs = 2.0j;o!g and logr^ = 2.3lo!g (r^ in light days) 

for a = 0.05 and a — 0.1 respectively. Taking into ac- 

1/2 

count that oc (Af) ' where (M) is the mean mass 
of the microlenses, our results for the mean size of the 
accretion disk of quasars are = 7.4^4 i((M)/Mq)^/^ 

light days for a = 0.05 and r,. = 10tll{{M) /Mq)^'^ 
for a = 0.1. When scaled for a mean microlens mass 
of (M) = O.3M0, representative of the stellar popula- 
tions in galaxies, the sizes of the accretion disks become 
rs — 4.012.2 light days and = 5.5^3;^ light days for 
a = 0.05 and a = 0.1, respectively, at our average wave- 
length (A)=1736A. The intrinsic dispersion in A in be- 
tween the sources implies an additional dispersion in Ts 
of about 20%. When taking this effect into account in the 
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+2.4 
-3.1 



and 



= 5.5 



+3.1 
3.7 



error bars, we finally have r, 
light days for a = 0.05 and a = 0.1, respectively. 

Using equation (3) and a Bayesian approach with a lo- 
garithmic prior on (see e.g. Mediavilla et al. 2011a), 
we also made individual estimates of for each one of 
the 19 quasars in our sample. For those objects for which 
there are several image pairs, the probability distribu- 
tions for the different pairs were combined into a single 
probability function for that object. The resulting values 
for {M) — 0.3Afo and a = 0.05 are presented in Table 
[T] and in Figure 3. It is worth noting that the average 
value of these individual estimates {{rg) = 4.2 ± 4.4 light 
days) is very close to the one obtained with the joint 
likelihood analysis (although the use of a linear prior 
would increase the individual size estimates by 80%). 
Note that the average size shrinks noticeably (by ^^30%) 
if the pairs with \AmB\ < 0.1 are excluded from the 
sample. Thus we confirm that favouring the selection of 
quasars with noticeable microlensing introduces a clear 
bias towards smaller sizes (Kochanek 2004, Eigenbrod et 
al. 2008, Blackburne et al. 2011a). In Figure|3]we com- 
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TABLE 1 

Sizes for individual objects at 1736 A 



Object 


rs(1736)A 




Line 


HE0047-1756 


4.0+«-^ 


1.38 


(CIV) 


HE0435-1223 


4.4l|i 


0.50 


(CIV) 


HE0512-3329 








SDSS0806+2006 








SBS0909+532 


i.stl-l 


1.95 


(H^) 


SDSSJ0924+0219 


10.1-4.6^ 


0.11 


(Mgll) 


FBQ0951+2635 


1.8+2| 


0.89 


(Mgll) 


QSO0957+561 




2.01 


(CIV) 


SDSSJlOOl+5027 








SDSSJ1004+4112 




2.02 


(CIV) 


QSO1017-207 


q r, + 6:02 

'^■^-2.12 


1.68 


(CIV) 


HEl 104- 1805 


1 0+"-^ 

^■^-0.5 


2.37 


(CIV) 


PG1115+080 


9 0+2.6 


0.92 


(CIV) 


SDSSJ1206+432 








SDSSJ1353+1138 


7.5-4.4" 






HE1413+117 


9 7+2.7 
^- ' -1.3 


0.26 


(CIV) 


BJ1422+231 


20.017/2 


4.79 


(CIV) 


SBS1520+530 


2.4tli 


0.88 


(CIV) 


WFIJ0233-4723 


3.6tli 







Note. — Sizes in light days for individual objects at 1736 A 
obtained from microlonsing magnifications of image pairs. Black 
hole masses are taken from the suimnary in Mosquera & Kochanek 
(2011). Emission line used for the mass estimate is indicated in 
fourth column. 

* These values have been calculated for Ams = and are lower 
limits rather than true estimates 

pare our estimates with previous results. Estimates from 
the present vi^ork are indicated with filled symbols, and 
estimates from previous works with open symbols. We 
have taken black-hole masses from Table 1 in Mosquera 
& Kochanek (2011), which they take from Peng et al. 
(2006), Greene et al. (2010) and Asset ct al. (2011). 
The thin line corresponds to the theoretical prediction 
for the thin disk model at 1736 A (with L/Le = 1 and 
T] = 0.1). The thick line corresponds to the empirical 
fit by Morgan et al. (2010) scaled to (A)=1736 A using 
P = 4/3 

Our estimates are in good agreement with previous 
results (Blackburne et al. 2011a, Mcdiavilla et al 2011a, 
Mosquera et al. 2011, Munoz et al. 2011, Morgan et al. 
2010, Floyd et al. 2009, Anguita et al. 2008, Poindexter 
et al. 2008, Eigenbrod et al. 2008, Bate et al. 2008, and 
Kochanek 2004) and show a discrepancy with respect to 
the predictions of the thin disk model of approximately 
a factor of 5. 

4. CONCLUSIONS 

We performed a statistical analysis to estimate the av- 
erage size of the accretion disks of lensed quasars from 
the large sample (MED09) of microlensing magnification 
measurements. We find that most information on the 
size of the accretion disk is contained in the amplitude 
of the microlensing magnification and fairly independent 
of the amount of chromaticity. 

From a statistical analysis using measured microlens- 
ing magnification strengths from 27 image pairs from 19 
lensed quasars in the sample of MED09 we measured the 
average size of the accretion disk of lensed quasars at a 
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Fig. 3. — Size of the accretion disk at 1736 A as a function of 
the mass of the black hole (from Mosquera & Kochanek 2011) for 
{M) = 0.3Mq . Results from our joint analysis arc the solid trian- 
gle and square for a=0.05 and a=0.1 respectively Our estimates 
for individual objcts arc shown as filled diamonds (or upward ar- 
rows for lower limits). Open triangles arc from Blackburne et al. 
(2011a) for objects with available black hole masses in Mosquera 
& Kochanek (2011). The asterisk is the estimate for SBS0909 from 
Mcdiavilla et al. (2011a). Estimates for HE1104 by Muiioz et al. 
(2011) and Poindexter et al. (2008) are shown as an open square 
and an X respectively. The estimate for HE0435 by Mosquera et 
al. (2011) is shown as an open circle. Upper half circle, lower half 
circle and star show the estimates for QS02337 by Anguita et al. 
(2008), Kochanek (2004) and Eigenbrod et al. (2008) respectively. 
The estimate for Q0957 by Hainlinc et al. (2012) is shown as a 
bow tic. Downward arrows indicate the upper limits for MG0414 
by Bate ct al. 2008 (right) and for SDSSJ0924 by Floyd ct al. 
2009 (left). Thin line is the prediction from the thin disk model 
(with L/Le = 1 and r)=0.1). The thick solid line is the fit from 
Morgan et al. (2010) with their measurements shown as pluses. 
The dashed lines showing the limits due to the uncertainty in the 
y-intercept. 

rest wavelength of (A) =1736 A. For a typical mean mass 
in the stellar population of the lens of (M) = O.3M0, and 
taking into account the intrinsic dispersion in A, we find 
rs = 4:.0tl[t ^™d rs = 5.5+^7 light days for a = 0.0-5 and 
a = OA respectively. These estimates are in agreement 
with other studies (Blackburne et al. 2011a, Mcdiavilla 
et al 2011a, Mosquera et al. 2011, Munoz et al. 2011, 
Morgan et al. 2010, Floyd et al. 2009, Anguita et al. 
2008, Poindexter et al. 2008, Eigenbrod et al. 2008, 
Bate ct al. 2008, and Kochanek 2004) and again find 
that disks arc larger than predicted by the thin disk the- 
ory. We have also estimated the sizes of the 19 individual 
objects in the sample to find that they are in very good 
statistical agreement with the result from the maximum 
likelihood analysis. 
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